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In this letter a generalization of the BEC-BCS crossover theory to a multicomponent superfluid is 
presented by studying a three-species mixture of Fermi gas across two Feshbach resonances. At the 
BEC side of resonances, two kinds of molecules are stable which gives rise to a two-component Bose 
condensate. This two-component superfluid state can be experimentally identified from the radio- 
frequency spectroscopy, density profile and short noise measurements. As approaching the BCS 
side of resonances, the superfluidity will break down at some point and yield a first-order quantum 
phase transition to normal state, due to the mismatch of three Fermi surfaces. Phase separation 
instability will occur around the critical regime. 



Introduction. During recent years, one major progress 
in quantum gases is achieving the crossover from a Bose 
condensation of molecules to a fermionic BCS superfluid 
in two-species fermionic gases across Feshbach resonance. 
Meanwhile, the multi-component Bose gases also attract 
considerable interests, mainly because of the effects re- 
lated to internal phase coherence. Here we arise an inter- 
esting question that how the crossover process happens if 
the molecular BEC is a multi-component one, and it will 
be intriguing to achieve a fermionic multi-component su- 
perfluid which should be more attractive than its bosonic 
counterpart. 

Lucky, this question is not only an academic one, but 
also close relates to current experiments on lithium and 
potassium gases. Recently an accurate measurement of 
the scattering lengths between the lowest three hyperfine 
spin states of lithium, denoted as |2) and |3) respec- 
tively, finds that there is a Feshbach resonances between 
|2) — |3) at Sl.lmT followed by another resonance be- 
tween |1) - |2) at 83.4mT, as shown in FigHQ. Similar 
feature has also been found between three hyperfine spin 
states of potassium, where two resonances are located at 
20.21mT and 22.42mr respectively 0. Note that both 
the |1) — |2) molecule and the |2) — |3)molecule are sta- 
ble at the BEC sides of two resonances pj, we consider the 
case that the numbers of atoms in different species satisfy 
N2 = N1+N3, therefore we have a two-component molec- 
ular condensate deeply in the BEC side. The main pur- 
pose of this letter is to present a mean-field theory to de- 
scribe two-component BEC-BCS crossover using lithium 
gas as an example 0- 

On the other hand, it is also interesting to look at the 
system from the BCS side. Note that the normal state 
has three different Fermi surfaces, how the system copes 
with gaining paring energy and Fermi surface mismatch 
is a long-standing problem for both condensed matter 
and high energy physics. Recent experiments on two- 
species mixture with population imbalancej^ have re- 
vealed many interesting phenomena such as superfluid- 
normal transition and phase separation, and also caused 
a considerable amount of theoretical interests j^. An- 
other point of this letter is to demonstrate that similar 



effects will also occur in the three-species mixture as it 
approaches the BCS side. A global phase diagram con- 
structed by the mean-field theory is shown in Fig^ 

Before proceeding, we would like to remark on the sta- 
bility of the system against atom loss, which is the major 
concern of experiments. There are three sources of atom 
loss: atom-molecule inelastic collision, molecule-molecule 
inelastic collision and the Efimov states. In order to re- 
duce the atom-molecule inelastic collision, one can ini- 
tially prepare the system deeply in the BEC side where all 
the atoms are bound and the number of unbound atoms 
is very few, and then adiabatically tune to the resonance 
regime. The inelastic collision between molecules could 
be largely suppressed by the Pauli exclusion principle 
as all molecules contain species |2). As for the Efimov 
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FIG. 1: An illustration of the scattering lengths ai2 and a23 
as a function of magnetic field and the phase diagram for the 
three-species mixture with A^2 = A''i -I- N^. The relation be- 
tween scattering lengths and magnetic field is obtained from 
the measurement of Ref.Q. In this range of magnetic field ai3 
is much smaller other two, and is not shown here. Restricted 
in homogeneous states, a first order quantum phase transition 
form the two-component superfluid state to the normal state 
is found above the first resonance. Incorporating imhomoge- 
neous states, the phase separated state will take over around 
the critical magnetic field. The locations of phase boundary 
will depend on the density and the concentration. 
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states, it is unlikely that it will significantly affect the 
three-body loss rate in this case, (see detailed discussion 
in footnote. and Ref.|3). Although the stability of this 
system remains to be seen in future experiments, qual- 
itative speaking it is very promising that it can have a 
reasonable long life time for experimental study, and we 
hope the theoretical discussed below will motivate more 
research in this system. 

Mean- field Theory of Two- component Super fluid. Fol- 
lowing the idea of crossover theory [^, we write down a 
BCS-like wave function which contains pairing between 
both |1) - |2) and |2) - |3) 10] 

l*> = n (^k + ^^kV-LV-lki + ^k^LV^ks) |0>- (1) 

k 

Note that the interaction between |1) — |3) is much weak 
in the range of magnetic field from 78mT to 88mT, com- 
paring to those between |1) — |2) and |2) — |3)P,|3], we 
neglect the pairing between 1 1) — 13) because both of them 
try to pair with |2) first for energetic consideration, and 
also because the critical temperature to achieve |1) — |3) 
pairing is too low to reach in current experiments. 

Here we would like to emphasize the symmetry and 
order parameters of this state. First, there are two in- 
dependent pairing order parameters, Ai = {''P^^A) ^'^'i 
^2 = (V'^V's)- Furthermore, as both species |1) and |3) 
pair with species |2), the operator tpsi^ii which in fact 
converts a |2) — |1) pair into a |2) — |3) pair, acts as a 
Josephson tunneling between two components, and the 
wave function Eq.(^ automatically gives another order 
parameter rj = {ip^ipi) which is related to the relative 
phase between two components. 

The Hamiltonian under consideration is written as 

W = X! ^kaV'LV'ka + X! 52iV'k2V'lkiV'-k'iV'k'2, (2) 
k(7 k,k',i 

where eko- = h?k'^ /{2m) — /i^, cr = 1, 2, 3 and i = 1, 3. As 
in the two-species mixture, g2i is related to the scattering 
lengthes via 1/521 = m/(A'Kh^a2i) — X]k™/(^^^^)- ^^^^ 
interaction between |1) — |3), which is neglected here, can 
be turned on as perturbation in a more detailed study 
elsewhere, and it will not affect the qualitative features 
discussed here. 

With the Hamiltonian Eq.® the free-energy for the 
quantum state of Eq. is given by 

^= E(eka|l^k|'+ekbh«k|') 
k 

+ X! "k^k^k' "f^k' + 523WkWkUk' If'k' , (3) 
kk' 

where ^ka = eki + ek2 and ^kb = ek2 + Cks- The con- 
straint |ukp -I- \v\s\^ + |wkP = 1 can be imposed by a 
Langrange multiplier X]k^k(|ukP + |wkP + |wkP - 1). 
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FIG. 2: The solution to the mean-field equations. (A-B): As 
the change of magnetic field, the evolvement of two pairing 
order parameters Ai (solid line) and A2 (dashed line) (Fig. A), 
and the chemical potential /ia = /ii + (solid line) and fib = 
M2 +M3 (dashed line) (Fig.B). (C-D): The momentum distribu- 
tions for 1 1) (dashed line), |2) (solid line) and |3)(dashed-dotted 
line) at 80.5mr(Fig.C) and 84mr(Fig.D). Here we have cho- 
sen n2 = 1.14 X 10" cm"^ and : A^i : A^3 = 1 : 0.5 : 0.5. 
The two arrows in (A) and (B) indicate the locations of Fesh- 
bach resonances, fcp is the Fermi momentum for the majority 
species |2), and /io ~ h^kp/{2m). 



Minimizing the free-energy with respect to Wk, fk and 
Wk gives dT /du\s_ — dT ldv\^ = dT /dw\^ = 0, which are 

Ak Ai A2 \ / "k \ 

AJ Ak - Cka ^'k = 0, (4) 

A^ Ak - Ckb J \Wk J 

where Ai = -g2iJ2k'^kV^ and A2 = -323 Z^k "k^i^- 
Thus Ak satisfies the equation 

Ak - ^kA^ + BkAk + Ck - 0, (5) 

where ^k = Cka + ^kfc, Sk = ^kaCkh - |Aip - IA2P and 
Ck = |Ai|2^kb + |A2pCka- The lowest solution of Eq.© 
is 

Ak - ^{Ak - 2^JaI - 3BkCos[(7r - 0k)/3]}, (6) 

where 61, = arctan[3V3i?k/(2A^ - 9AkBk - 27Ck)l and 
Kk = AlBl - ABl + 4^3 Ck - IS^lkSkCk + 27^2 11]. 

For the wave function satisfying Eq. I^J , the free-energy 
is given hy T = Y,k^^~ l^iP/521 - 1A212/523- H can 
be verified that the long wave length behavior of Ak is 
precisely — (JAij^ + \/S.2\'^)m/ {h^k'^), therefore the diver- 
gency in the summation over Ak can be exactly canceled 
out by the renormalization terms in 321 and (723 • Fur- 
thermore, dJ- /d/S* =0 [j — 1, 2) yield two coupled self- 
consistency equations. Differential Eq.® with respect 
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to Aj, one finds that 9Ak/5A* satisfies 



(3X1 - 2AkAk + Bk) 



5A* 



Thus, dT/dAl gives 
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Ak-6 



kb 



3A2 - 2AkAk + Bk 



Ai(Ak-ek;,). (7) 



(8) 



and similarly we have 
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3A^ 



2AkAk + Sk h?k^/m 



(9) 



Apart from increasing of gaps, another feature of the 
BEC-BCS crossover is the shift of chemical potentials, 
which is determined with the help of number equations. 
From Eq.Q) we can get Uk/wk — (■fka — Ak)/A^ and 
u^kK = [Ak(Ak - ^ka) - |Ai|2]/(AtA2). Therefore, 
two number equations are given by Ni = |wkp = 
Ek(l/(1 + l^kMP + kkMP) and iV3 = Ek l^kP = 
Ek(kk/ukP/(l + luk/wkP + kk/wkP) respectively. 

We solve these four coupled equations numerically, and 
the results are shown in Fig 13 which exhibit the com- 
mon features of BEC-BCS crossover. From Fig[2lA-B), 
one can see that the gaps increase and the chemical po- 
tentials decrease as the magnetic field changes from the 
BCS side to BEC side. Fig|2i;C-D) displays the momen- 
tum distributions for three different species. One can see 
that the momentum distributions become much broader 
at the BEC side comparing to those at the BCS side, 
which indicates that the real space sizes of two kinds of 
Cooper pairs are much smaller. This solution also allows 
one to calculate the energy of this superfluid state E^f- 

Quantum Phase Transition for Homogeneous System. 
As we have mentioned, in the normal state the Fermi 
surface of species |2) is larger than those of species |1) 
and 1 3) for the concentration N2 = Ni + N3, it will cost 
a lot of kinetic energy of two minority species to form 
pairs. Thus one expects that there is a quantum phase 
transition form this two-component superfluid state to 
the normal state as approaching BCS side. 

At mean-field level, we simply take the normal state 
as non-interacting normal state, and its energy i?N is the 
total kinetic energy of three species. In Fig|3{A-B) we 
plot the energy difference SE = i?N — Esf for different 
densities and concentrations. A quantum phase transi- 
tion occurs when 6E becomes negative, i.e. the normal 
state has lower energy. The critical point is usually lo- 
cated around the first Feshbach resonance, and its exact 
location depends on density and concentration. We find 
that the critical point will be pushed toward the BCS 
side either when the densit y o f atoms increases or when 
the radio A^i : increases [l^- We also notice that one 
can not connect the quantum state of the form Eq.QJ 
to a normal state by continuous varying the parameters. 
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FIG. 3: (A-B) SE as a function of magnetic field. (A): N2 : 
iVi : iVs = 1 : 0.5 : 0.5 and kpao = 1 x 10"'' (solid line), 
fcpao — 0.5 X 10"'' (dashed line), fcpao — 2 x 10"'' (dashed- 
dotted line). (B) kpao = 1 x 10"* and iV2 : iVi : iVs = 1 : 0.5 : 
0.5(solid line), iV2 : iVi : iVa = 1 : 0.4 : 0.6(dashed line), N2 : 
iVi : iVa = 1 : 0.6 : 0.4(dashed-dotted line). (C) Comparison 
of the energy between the phase separated state(dashed line), 
the two-component superfluid state(solid line) and the normal 
state(dotted line), fcpao = 1 X 10"" and N2 : Ni : N3 = 
1 : 0.5 : 0.5. For all figures the unit of energy density is 
fp = kp / {20n^ m) , fcp is the Fermi momentum of the major 
species |2) in its normal state, and ao — 0.0529177nm. Two 
arrows indicate the locations of Feshbach resonances. 



therefore this transition should be a first order quantum 
phase transition. 

Phase Separation Instability. There are possible phase 
separation instability in both normal state and the 
two-component superfluid state. Approaching the first 
resonance from the normal state, ai2 becomes larger 
while a23 is still not larger, the energy may be low- 
ered if the normal state separates into a phase sepa- 
rated state constructed as follows: in one region it is 
nk("k + ^'kV'LV'lki) ripLo V'pslO) and in another region 
it is nk=o V'k2 ripLo V'pslO). The energy of this state is 
shown in Fig|3JC) for a typical density and concentra- 
tion. We find that around the critical magnetic field the 
energy of phase separated state is lower than both two 
homogeneous states. 

In the superfluid phase k = d^T/{d\Ai\'^d\A2\^) is al- 
ways positive, indicating the interaction to be repulsive. 
Hence, one need to consider whether the repulsion is so 
strong that lead the two-component condensate spatially 
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separating into two condensates, with |2) — |1) molecules 
staying in one side and |2) — |3) molecules staying in an- 
other side. Here we also calculate the energy of this type 
phase separated state, and find that the energy is very 
close to, but usually sHghtly higher than, the energy of 
the homogenous state in all range of magnetic field. 

Combining the discussion on the homogenous states 
and the phase separation instability, we draw the con- 
clusion on the phase diagram illustrated in Fig^ We 
remark that the effects missing in the mean-field the- 
ory, including the interaction energy of normal state and 
quantum fluctuations in superfluid phase, will give quan- 
titative corrections to this phase diagram. However, the 
qualitative feature, which is guaranteed by the physical 
understanding at two ends, should holds. 

Experimental Signatures of the Two-component Super- 
fluid. At the end of this paper, we point out some ex- 
perimental signatures to reveal the unique features of the 
two-component supefluid state. 

To detect the coherence between two components, one 
can look at the radio-frequency spectroscopy, which can 
be described by iJrf = V'sV'i + h.c. and has already been 
widely used in fermion experiments^3j| . Applying r-f field 
in this superfluid state will induce a time-perodic oscilla- 
tion of the atom number in both species |1) and |3), which 
can be a direct evidence of phase coherence between two 
components. 

A direct signal to distinguish this state from the nor- 
mal state is an in situ measurement of the density profile. 
In the superfluid state the density profiles of different 
species have to satisfy n2{r) — ni{r) -|-n3(r) everywhere, 
and the Thomas-Fermi radii are the same for all species. 
In contrast, in the normal state the Thomas- Fermi radius 
of the majority species is larger than those of two minor- 
ity species, and n2(r) does not equal to ni(r) -|-n3(r) ev- 
erywhere. For non-interacting normal state, at the center 
of the trap, n2/(ni + ^3) = N^^^/{N^^^ N^^^) < 1. 

Analyzing the noise correlation in time-of-flight im- 
age allows one to measure the second-order correlation 
function GafjiKk') = (nQk«/3k') - ("-ak) ("/3k')'l^- In 
this state e2i(k,k') = |ukVkP^(k k') and e23(k,k') = 
|wkWkP<5(k -I- k'), which shows strong pairing correlation 
between opposite momentums. Here, one remarkable fea- 
ture is the effective exclusive correlation between |1) and 
|3), and it shows up in tJi3(k, k') ~ — |wkWkP^(k — k') 
as a dip at equal momentum. This arises from the fact 
that the atom in the quantum state |2, — k) can only pair 
with either |f,k) or |3,k), thus the quantum state with 
momentum k can not be occupied by both |1) and |3). 
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